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Abstract 

The master field for a subclass of planar diagrams, so called rainbow diagrams, 
for higher dimensional large N theories is considered. An explicit representation 
for the master field in terms of noncommutative random variables in the modified 
interaction representation in the Boltzmannian Fock space is given. A natural inter- 
action in the Boltzmannian Fock space is formulated by means of a rational function 
of the interaction Lagrangian instead of the ordinary exponential function in the 
standard Fock space. 
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1 Introduction 



The problem of summation of all planar diagrams in higher dimensional space-time is 
still out of reach. Its solution is closely related with problem of finding the leading 
asymptotics in matrix models for large N and may have important applications to the 
hadron dynamics [pi 0, 0]. Summation of planar diagrams has been performed only in 
low dimensional space-time Q. 

One can write a closed system of equations for invariant correlation functions in the 
large N limit, so called the planar Schwinger-Dyson equations, for arbitrary dimension of 
space-time. In the early 80-s it was suggested || that there exists the master field $(x) 
such that the correlation functions for this field are equal to the large N limit of 

invariant correlation functions for matrix models, 

jirn < tr(M(x n )...M( Xl )) >=< 0|$(x n )...$(x 1 )|0 > (1.1) 

It was suggested that satisfies the equation |J 

[i^ + m{x)]\o>=o (i.2) 

where S is an action and II and $ are the subjects of the relation || 

[U(x), <%)] = i5 {D \x - y)|0 >< 0| (1.3) 

An operator realization of this algebra proposed in || has used the knowledge of all 
correlation functions. This was considered as an evident drawback of such approach. In 
that time it was also discussed the problem of finding a generating functional reproducing 
the planar Schwinger-Dyson equations or equations fli~2]), (TTTH) . It was pointed out that 
the generating functional cannot depend on one commutative source 0. It was proposed 
to use some auxiliary fermionic fields to reproduce the planar Schwinger-Dyson equations 
||. For gauge theories the suitable generating functional is nothing but a functional 
on paths, i.e. the Wilson loops, that satisfies the Makeenko-Migdal equation ||. The 



stochastic equation for large N master fields was proposed in |"U 



Recently it has been a reveal of an interest to the problem of constructing the master 



field for planar graphs. One of origins for this are mathematical works [|11|, E| [13| devoted 



to non- commutative probability, for a review see |1J, Singer has advocated that the 
essential difficulty of the large number of degrees of freedom in higher dimensional large 
N matrix models is dealt with finding master fields which live in " large" operator algebras 
such as the type ll\ factor associated with free group. In the recent paper by Gopakumar 
and Gross [|16| the basic concepts of non-commutative probability have been reviewed and 
applied to the large N limit of matrix models. They stress that if one can solve a matrix 
model then one can write an explicit expression for the master field as an operator in a well 
defined Hilbert space. Douglas also proposed to use ideas of non-commutative probability 
to the large N stochastic approach [jrTJ. The explicit construction of the master fields for 



several low- dimensional models including QCD2 has been given [ |16| 

However an effective (i.e. without the knowledge of correlation fonctions) operator 
realization for the master field for all planar diagrams in higher dimensional space-time 
is still unknown. Therefore it is worth to try to find an effective operator realization of 
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Figure 1: Wick theorem for rainbow diagrams 



the master field for some subset of the planar diagrams. The goal of this letter is to 
construct an explicit operator realization of the master field for rainbow graphs. Rainbow 
graphs form a subset of planar graphs. It is turn out that rainbow correlation functions 
may be obtained by average of the fields with Boltzman statistics. The construction 
does not require as input correlation functions. More exactly we show that to get a 
closed set of equations for correlations functions for a model with an interaction in the 
Boltzmannian Fock space one has to deal with a modified interaction representation. 
This new interaction representation involves not the ordinary exponential function of the 
interaction but a rational function and will be given by the formula 

< <t>{x m )..4{xi)- 1 > (1.4) 

1 - 9 1 dyV int {(j){y)) 

The paper is organized as follows. In Sect. 2 we present the Schwinger-Dyson equations 
for rainbow diagrams. Sect. 3 contains a necessary information about the Boltzmannian 
Fock space. We argue also that to get a closed set of the Schwinger-Dyson type of equation 
we have to deal with the interaction representation in the Boltzmannian Fock space in the 
form ( |1.4|) . We also show that the corresponding Schwinger-Dyson equations reproduce 
the Schwinger-Dyson equations for rainbow diagrams. 

2 Rainbow Diagrams 

Let us consider the correlation functions of the form 

< V(x n , ...x x ) >= < tr ( M(x n )...M( Xl )) > (2.1) 

< 0(M) >= - f 0(M)exp{-S[M}}dM, (2.2) 
S[M] = J dx[^ tr ( M(- A +m 2 )M) + ^^ M \ x )} ( 2 -3) 



< • > means 



where 



Here M(x) is N x N matrix function. We assume all necessary regularizations. For 
our purpose it is essential that the regularization is such that the free propagator is 

< M ij (x)M fi/ (y) ><°>= 5 iif 5tfD(x - y), (2.4) 
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Figure 2: Schwinger- Dyson equation for planar graphs in the large N limit 

(- A+m 2 ) x D(x -y) = 5 {D \x -y). (2.5) 

We shall consider the external lines corresponding to global invariant Green func- 
tions as the lines corresponding to generalized vertex. The rainbow diagrams for 
< V(x n , ...x\) > are defined as a part of planar non- vacuum diagrams which are topologi- 
cally equivalent to the graphs with all vertexes lying on some stright line on the the right 
of generalized vertex and all propagators lying in the half plane. The rainbow diagramms 
are illustraited for the M 3 - interaction on Fig.|l|, where vertexies are drawn by solid double 
lines and all contraction (propagators) by double dash lines; in that follows we will use 
also solid lines for propagators. 

To write down the rainbow Schwinger- Dyson equations one has to select rainbow 
diagrams from the both hand sides of the planar Schwinger-Dyson equations. In the large 
N limit the planar Schwinger-Dyson equations have the form 

(- A +m 2 ) Xl G n (x n , ...xi) = gG n+2 (x n , ...x t+1 ,x h x h x h x\- X , —X\) (2-6) 

+ 5(xi — Xj)G;_j_i(Xj_i, ...Xi + i)G n+i -i^i(x n , ...Xt+i,Xi-i, ...Xi) 

i<l 

+ ^25(xi - ...x i+ i)G n+ i-i-i(x ri , ...x i+ i,xi- 1} ...Xi), 

Ki 

where 

G n (x n ,... Xl ) = Jiin < tx{M{x n )...M{x x )) > (2.7) 

The planar Schwinger-Dyson equations (|2.6| ) are written for the case of quartic interaction 
( |2.3| ) and they are symbolically presented on Fig.|2]. Now let us consider a modification 
of the right hand side of (|2.6|) for correlation functions corresponding to the rainbow 
diagrams 

i im ^nT^ < tr ( M(x n )...M( Xl )) > rb = W n+2 (x n+2 , ...x x ) (2.8) 
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Figure 3: One term in the Schwinger- Dyson equation for rainbow graphs 

There are modifications in the term representing the interaction and also in the Schwinger 
terms. Indeed, let us consider all possible contractions of a given point x m with a vertex 
v of the rainbow diagrams. We have to distinguish the cases of odd and even m. For the 
even m = 21 on the left of this vertex v (see Fig.[3]a) one has subgraphs corresponding to 
rainbow diagrams of the correlation function (Fig||b) 

< (M(x 2l )M(x 2l _ 1 )...M(x 1 )(V mt ) k %, > (2.9) 

or 

< (M 3 (x 2l )M(x 2l _ 1 )...M(x 1 )(V mt ) k >) n > (2.10) 

By using 

r 

< (M(xi)...M(x k ))jji >= < tr ( M( Xi )...M(x k ) > (2.11) 

these terms reproduce W 2 i{x 2 h x 2 ;_i...Xi) and W 2l 1 +2 (x 2 i,x 2 i,x 2 i,x 2 i^i...Xi), respectively. 
The rest of the diagram Fig.|3|a corresponds to W k ^_ 2l+2 (x nj ...x 2 i+i,x 2 i, x 2 i), k 2 = k—k\ — l. 

One has also to modify contributions from the Schwinger terms, since only one corre- 
lator (in our case the correlator corresponding to G n -i+i-i) can contain the interaction. 
Finally we get the following system of equations 

(- A +m 2 ) x . 2l W n (x n , ...xi) = giyV^x-n.x-n-x, ...xi)W n ^ 2 i +2 (x n , ...x 2l+1 ,x 2i ,x 2 i) 
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+ W 2l+2 (x 2h %2h %2U %2l-l, ...Xi)W n _2l{x w ...X 2 l+l)) 
+ $( x 2l ~ Xi)W n ^ 2 i +l - 3 (x n , ...X 2 l+1, Xi-i, ...Xi)W 2 il i+l (x 2 l-l, ...Xi+x) 



i<2l 

r(0) 



+ ^2S(x 2 i - Xi)W n+2 i-i-3(x n , ...x i+ll x 2 i-i, ...xi)W;_ 2l+1 (x i+1 , ...X21+1); (2.12) 



2l<i 
2 



(- A+m ) X2l+l W n (x n , ...xi) = g(W 2 i(x 2 i,X2i, ■■■Xi)W n - 2 i+2(x n , ...X21++2, x 2 i+x, x 2 i+i, x 2 i+ 



+ W 2 l + 2(X21+1, X21+1, X 2 l, ...Xl)W n -%{x m ...X21+2X21-1)) 

+ S ( X 2l ~ Xi)W n ^2l+i-3(x n , —X21+I, Si-l, ...X 1 )W 2 iL i+1 {x2l-l, ...X i+ i) 
i<2l 

+ J2 <K X 2Z ~ Xi)W n+2 i-i^ 3 (x n , ...x i+1 ,x 2 i-i,x 1 )w£% l+1 (x i+ i, ...x v+ i) (2.13) 

Note that for planar correlation functions it is enough to write down the equation 
differentiated only on x\ since the function G(x n , ...x\) is invariant under cyclic permuta- 
tion. The rainbow correlation functions loss this property and we have to write differential 
equations for all points x\. In the the right hand side of ( 2.12|) and (|2.13| ) enters t he free 
rainbow correlation functions W®(x n , ...X\). They satisfy the system of differential equa- 
tions being the subject of equations ( |2.6|) with g = 0. In distinguish of the interaction case 
the free rainbow correlators posses the property of invariance under cyclic permutation. 
This is due to the fact that free planar correlators coincide with free rainbow correlators, 
i.e. = To find it is enough to apply the Wick theorem as it is shown on 

Fig.|a. 



3 The Master Field for Rainbow Diagrams 
3.1 Free Field Theory 

Let us consider an algebra generated by operators A(p) and A + (p) satisfying the relations 

A(p)A + {q) = ^(p-q). (3.1) 

One can realized this algebra in a space which is an analogue of the usual Fock space 
[Ll|. This space is generated by the vacuum vector |0 >, A(p)|0 >= and n-particle 
states of n non-identical particles, 

\p u ... Pn >= A+( Pl )...A+(p n )\0 > (3.2) 

There is no symmetization or antisymmetrization as in the Bose or Fermi cases. We shall 
call this Fock space the Boltzmannian Fock space (it is also called the free Fock space). 
One defines 



!.r) „+(,-) + n-u) = / J^ (A+(p)t*° + A(p)e*") (3.3) 



and therefore 



< om<t>(y)\Q >= D ^-V) = j^y5j ^T^^ {X - y) (3-4) 
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Figure 4: a) Wick theorem for rainbow graphs in free matrix theory; b) Wick theorem in 
the Bolzmannian Fock space 

To calculate the n-point correlation function one has to apply a Boltzmannian Fock space 
analog of the ordinary Wick theorem. The specific feature of the Wick theorem in this 
case is that for a given diagram one has not additional symmetry factors related with 
that an annihilation operator can be contracted with any creation operator on the right. 
In the Boltzmannian Fock space an annihilation operator can been contracted only with 
a nearest creation operator on the right. Therefore one sees immediately from the Fig^b 
that the correlation function 

< 0\ ( f)(x 2m )... ( f>(x 1 )\0 >= F(°2(x2m, -xx) (3.5) 

satisfies to the same equations as W 2m , and therefore fjm^ra; ...x\) X 2 rni ■ •■X\ ) , 

i.e. 

fim — L- < tr ( M(x 2m )..-M(xi)) ><°>=< 0\ ( f>{x 2m )... ( j>{x 1 )\0 > (3.6) 

N— >oo iV ' 

Let us make a few comments about an operator realization of the algebra (|1.3|) with 
7r satisfying the requirement 

<x)\0 >=\{~ A +m 2 ),0(x)|O> . (3.7) 

First of all note that the operator algebra (|1.3|) is not an unique algebra which follows 
from the free planar Schwinger- Dyson equation. Indeed, one gets the same equation from 
the operator relation 

(- A +m 2 ) x (j)(x) = ir(x), [vr(x), <f>{y)] = -i5 D {x - y)|0 >< 0| + K{x, y) (3.8) 

with K(x, y) being the subject of relations 

< 0\K(y, Xl )<P{x2)...<p{x n )\0 > + < Q\(f)(x 1 )K(y,x 2 )...4>(x n )\0 > +... (3.9) 

+ < 0\<j)(x l )...<j)(x n „ 1 )K(y,x n )\0 >= 
The simplest solution of ( J3.8[ ) and fl3.9p is given by 

7T(y) = % -{- A +m 2 ) y [<p + (y)\0 >< 0| - |0 >< 0\<f>-(y)], (3.10) 
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K{y,x) = A +m 2 ) y [-<f> + {x) ( j> + {y)\0 >< 0| - |0 >< 0|<TW» (3.11) 

+<p + (y)\0 >< 0\<p~(x) + <p + (x)\Q >< 0\<p-(y)] 

This gives a hint to write a following operator realization of the commutation relations 
( |1.3|) with 7r satisfying the requirement ( |3.7|) 

n(y) = A +m 2 ) y {<p + (y)\0 >< 0| - |0 >< 0\<f>-(y)+ (3.12) 
f] / rfzi(- A +m 2 ) Zl ... / efe n (- A +m 2 ) Zn [<P + (z 1 )..4 + (zn)<P + (y)\® >< 0|<T (^)...<T (^ 

n=l ^ ^ 

-0 + (^ 1 )... ( /) + (^)|O >< O|0-(y)0 + (^ n )...0 + (z 1 )]} 

For completeness let us present a known solution of equation ( |2.(j| ) for D = and 
(7 = 0, i.e equations 

TO— 1 

($ 2m ) = V ($ 2Z )($ 2m - 2 ^ 2 ) (3.13) 
z=o 

where $ = a + + a, aa + = 1. Let us denote c n = ($ 2n ), c = 1. Consider the generating 
function 

Z{g) = c + c l5 + ... + c n g n + ... = ( — — ) (3.14) 

One has 

Z(#) 2 = Cq + (c Ci + CiCo)^ + ... + (c C n + ... + CnCo)^™ + ... 

From (|3.13| ) one gets 

Z{gf = Cl + c 2 £ + ... + c n+ i# n + ... 
Therefore Z(g) satisfies the equation 

gZ{g) 2 = Z(g)-l. 

One has to take the following solution of this equation 

Z(g) = l^PEEL = i + C \g + ... + ^—C- n g n + 
2g n + 1 



from which we get 



c -< $ 2 « >= __ a 



1 ™ 2n! 



n + 1 2n n!(n + 1)! 



3.2 Interacting Theory 

We want to derive the Schwinger-Dyson equations for theory with interaction in the 
Boltzmannian Fock space. To find the form of interaction let us consider the following 
correlation functions 

Fg\x m , ... Xl ) =< 0\(f>{x m )... ( f>{x 1 )( [ dy x : {4>{yi)) A :)...( / dy k : (0(y fc )) 4 :)|0 >' (3.15) 



where <j>{x) is the free field ( p.3|) and / means that we do not take into account the diagrams 
with vacuum subgraphs. We draw all operators (f) on the staight line. The operators 4>(xi) 
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Figure 5: Some diagrams contributed to the Schwinger- Dyson equations in the Boltz- 
mannian Fock space 

corresponding to the external lines are represented by the circles and the operators <f)(yi) 
corresponding to the interaction vertices are represented by the filled circles on Fig.|5]. 

On Fig.^] we draw all possible contractions of a given external line with a given vertex. 
As we have mentioned above there is no here additional factors related with symmetry 
of graphs. Therefore one has not here the standard factor l/k\ in the k-th order of 
perturbation theory. This remark leads to an important observation that to get a set of 
equations for correlation functions in the Boltzmannian Fock space we have to consider 
instead of the usual exponential factor exp{Vi„ f } the rational function 

{1-V int }~\ (3.16) 

(compare with Z(g) ( |3.14| )). Therefore we introduce the following correlation functions 

F m (x m ,... Xl ) = Y,9 k F£\x m ,...x 1 ) =< U\4>{x m )..4{xi) * |0 >' (3.17) 

fc=0 1 Vint 

OnFig.!^^^/^^)) 4 . 

Examining all possible contractions of the point x 2 i we see that the first two graphs 
on Fig.|5| reproduce the first two sums on the right hand side of fl2.12|) . In the similar way 



one sees that for odd point the graphs reproduce the first two sum in the the right hand 
side of ( gTg ). 



Performing the normal ordering in the expression (J)(x m )...(f)(xi) we get contributions 
corresponding to the Schwingers terms in the correlation function (— A-\-m 2 ) X2l F m (x m , ...x\). 
This consideration proves that F m satisfies to the following equations 

(- A +m 2 ) X2l F n (x n , ...X!) = g{F 2 i(x 2 i, x 2 i-i, ...xi)F n _ 2 i+ 2 (a; n , ...x 2 i+i,x 2 i, x 2i ) 
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+F 2 l + 2{X2U x 2h %2h %21-Xi ■■■Xl)F n _2l(x n , ...X2I+1)) 
+ K X 2l ~ Xi)Fn-2l+i-3{Xn, -X21+1, ■ ■ ) -f^f-i+l { x 2l-l , •••^i+l) 



i<2l 

<0) 



+ X! ^( X 2« — Xi)F n+ 2i-i-3(x n , ...Xi+i, X21-1, ■■■Xi)F i _ 2l+1 (x i+ i, ...X21+1); (3.18) 

2l<i 

and the similar equations for the a^z+i- 

Comparing (|3.18|) with ( 2.12|) we see that F m satisfies to the equations for the rainbow 



diagrams. Therefore we get 

Ji'JSLw^ < tr(M(x m )...M(i 1 ))exp{ 3 | F J d„ tr ( M (y)) 1 } >J,- (3.19) 

<0 W x m )...^,) 1 _ g / w |0>^ 

Symbol < . > F denotes the vacuum expectation value in the ordinary Euclidean 
bosonic Fock space and < . > BF denotes nonvacuum diagramms in the Boltzmannian 
Fock space. 

In conclusion, a model of quantum field theory with interaction in the Boltzmannian 
Fock space has been considered. We have used the new interaction representation with a 
rational function of the interaction Lagrangian instead of the exponential function in the 
standard interaction representation. The Schwinger-Dyson equations were derived and it 
was shown that the perturbation expansion for the model corresponds to the summation 
of the rainbow diagramms. The quantum field with this interaction can be interpreted 
as the master field for the rainbow diagramms in the largs iV limit matrix model. The 
construction of the master field is effective in the sense that it is purely algebraic and 
doesn't require the knowledge of correlation functions of the theory with the interaction. 
Another aspects of quantum field theory in the Boltzmannian Fock space including the 
Minkowskian formulation are considered in 
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